#include <iostream>
#include <Eigen/Dense>

int main(){
    Eigen::Matrix3d t1_A, t1_B;
    Eigen::Array33d t1_a, t1_b;
    t1_A << 1,2,3,4,5,6,7,8,9;
    t1_B.setIdentity();
    // t1_B = Eigen::Matrix3d::Identity(3,3);
    t1_a << 1,2,3,4,5,6,7,8,9;
    t1_b = t1_B.array();
    std::cout << "t1_A = \n" << t1_A << "\n\n";
    std::cout << "t1_B = \n" << t1_B << "\n\n";
    std::cout << "t1_A * t1_B = \n" << t1_A * t1_B << "\n\n";
    std::cout << "t1_a = \n" << t1_a << "\n\n";
    std::cout << "t1_b = \n" << t1_b << "\n\n";
    std::cout << "t1_a * t1_b = \n" << t1_a * t1_b << "\n\n";

    return 0;
}